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Abstract 

Recently exact results for the complete fermionic two-loop contributions to the 
prediction for the W-boson mass from muon decay in the electroweak Standard Model 
have been published l s . This paper illustrates the techniques that have been applied 
for this calculation, in particular the renormalisation procedure and the treatment of 
IR-divergent QED contributions. Numerical results are presented in terms of simple 
parametrisation formulae and compared in detail with a previous result of an expansion 
up to next-to-leading order in the top-quark mass. An estimate of the remaining 
theoretical uncertainties of the M\y-prediction from unknown higher-order corrections 
is given. For the bosonic two-loop corrections a partial result is presented, yielding the 
Higgs-mass dependence of these contributions. 
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1 Introduction 




One of the most important quantities for testing the Standard Model (SM) or its extensions 
is the relation between the massive gauge boson masses, and Mz, in terms of the Fermi 
constant, G^, and the fine structure constant, a. This relation can be derived from muon 
decay, where the Fermi constant enters the muon lifetime, r M , via the expression 

Gl ml 

' - ;i) 

with F(x) = 1 — 8x — 12x 2 Inx + 8x 3 — x 4 . By convention, this defining equation is supple- 
mented with the QED corrections within the Fermi Model, Aq. Results for Aq have been 
available for a long time at the one-loop |2| and, more recently, at the two-loop level (3J. 
Commonly, tree- level W propagator effects giving rise to the (numerically insignificant) term 
3mj;/(5M{y) in eq. (JTJ) are also included in the definition of G^, although they are not part 
of the Fermi Model prediction. 

Comparing the prediction for the muon lifetime within the SM with eq. yields the 
relation 

where the radiative corrections are summarised in the quantity Ar jl]. This relation allows 
a prediction of M\y, to be tested against the experimental result for Mw- The current 
accuracy of the measurement of the W-boson mass, M^ p = 80.451 ± 0.033 GeV 0, will 
be further improved in the final LEP analysis and Tevatron Run II [BJ, each with an error 
of 5M W « 30 MeV. At the LHC, an error of 5M W « 15 MeV can be expected [7j, while a 
high-luminosity linear collider running in a low-energy mode at the W + W~ threshold could 
reach a reduction of the experimental error down to <5Mw ~ 6 MeV This offers the 
prospect for highly sensitive tests of the electroweak theory [9J, provided that the accuracy 
of the theoretical prediction matches the experimental precision. 

The quantum correction Ar has been under extensive theoretical study over the last 
two decades. The one-loop result jlj involves large fermionic contributions from the shift 
in the fine structure constant due to light fermions, Aq oc login f, and from the leading 
contribution to the p parameter, Ap, which is quadratically dependent on the top-quark 
mass m t , resulting from the top-bottom mass splitting [TO] . 

Ar^ = Aa - ^Ap + Ar rem (M H ), (3) 

with Syj = 1 — M^/M|. The remainder part Ar rem contains in particular the dependence 
on the Higgs-boson mass, Mr. 

Beyond the one-loop order, resummations of the leading one-loop contributions Aa and 
Ap have been derived They correctly take into account the terms of the form (Ap) 2 , 
(AaAp), and (AaAr rem ) at the two-loop level and (Aa) n to all orders. 
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Beyond the two-loop order, complete results for the pure fermion-loop corrections (i.e. 
contributions containing n fermion loops at n-loop order) are known up to four-loop or- 
der ^21- These results also include the contributions arising from resummation of Act and 
Ap. Recently, the leading three-loop contributions to the p parameter of 0(G^m®) and 
0{G 2 L a s ml) have been computed in the limit of vanishing Higgs boson mass ^3], but were 
found to have small impact on the prediction of the W mass. 

Higher order QCD corrections to Ar have been calculated at 0(aa s ) |T3] and for the 
top-bottom contributions at 0(aa 2 ) 15]. The 0(aa 2 ) contributions with light quarks in 
the loops can be derived from the formulae (29) - (31) in |TB] and turn out to be completely 
negligible. First results for the electroweak two-loop contributions have been obtained using 
asymptotic expansions for large Higgs |T7j and top-quark masses [T%lll9H2Uj . Concerning the 
expansion in m t , the formally leading term of 0(G 2 ^m^) [THIIT^] and the next-to-leading term 
of 0{G 2 l m 2 M^) j20] were found to be numerically significant for the prediction of the W mass. 
Since both contributions turned out to be of similar magnitude and of same sign, a more 
complete calculation of electroweak two-loop corrections to Ar without using expansions is 
desirable. 

As a first step in this direction, exact results have been obtained for the Higgs-mass 
dependence (e.g. the quantity M WiSubtr (M H ) = M W (M H ) - M W (M H = 65 GeV)) of the 
fermionic two- loop corrections to the precision observables [2T|. They were shown to agree 
well with the previous results of the top-quark mass expansion |2*2*j . 

For the bosonic two-loop corrections to Ar, the complete result is not available up to now. 
However, in Ref. (22] the effect of the bosonic terms up to O (a 2 ) on the relation between the 
MS and on-shell definition of the gauge boson masses has been studied. For this purpose the 
corresponding two-loop self-energies have been evaluated in the MS-scheme using large-mass 
expansions. 

This paper discusses the exact computation of all fermionic two-loop corrections to Ar 
which has been presented recently PQ. These include all two-loop diagrams contributing to 
the muon decay amplitude and containing at least one closed fermion loop (except the pure 
QED corrections already contained in the Fermi Model result, see eq. ([T])). Some typical 
examples are shown in Fig. ^ No expansion in the top-quark mass or the Higgs boson 
mass is made, so that the full dependence on m t and Mh as well as the complete light- 
fermion contributions at two-loop order are contained. Previously, corrections from light 
fermions have only been taken into account via resummations of the one-loop light-fermion 
contribution (the two-loop light-fermion contributions have been calculated within the MS- 
scheme in Ref. |24j ) . 

The result of pQ has been included in the Standard Model fits and the indirect deriva- 
tion of constraints on the Higgs boson mass performed by the LEP Electroweak Working 
Group [Sj. 

As a further step towards a complete two- loop result for Ar, a partial result is pre- 
sented for the purely bosonic electroweak two-loop corrections which yields the Higgs-mass 
dependence of these terms. 

The paper is organised as follows. Sections |2}B1 enlarge on the methods which were 
employed for the calculation of the fermionic two-loop corrections. While section 121 presents 
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Figure 1: Examples for types of fermionic two- loop diagrams contributing to muon decay. 



an overview over the techniques, the renormalisation procedure is explained in section El and 
the extraction of the QED corrections, which are already contained in the Fermi Model, 
is described in section 0] A discussion of the numerical results and remaining theoretical 
uncertainties due to unknown higher orders can be found in sections 03 and HI respectively. In 
section [7| the Higgs-mass dependence of the bosonic two-loop corrections is studied. Before 
concluding, an outlook to the situation for future colliders is given in section |H1 

2 Outline of the two-loop calculation 

This section presents an overview over the main features of the calculation. 

Since the definition of the Fermi coupling constant according to eq. contains QED 
corrections in the Fermi Model summarised in the quantity Aq, the corresponding contri- 
butions have to be identified and extracted from the Standard Model computation of muon 
decay in order to arrive at the quantity Ar. As shown in section 0] below, all IR-divergent 
loop contributions in our calculation are already contained in Aq, contributing to Ar are 
IR-finite. As a consequence, after extraction of the Fermi Model contributions, it is possible 
to neglect the masses and momenta of the external particles, thereby reducing the generic 
diagrams contributing to the muon-decay amplitude to vacuum diagrams. 

For the renormalisation the on-shell scheme is used throughout. It entails, in addition, 
the evaluation of two-loop two-point functions with non-zero external momentum, which is 
technically more involved. However, it should be noted that the evaluation of this type of 
integrals is generally necessary in all renormalisation schemes if the result shall be related 
to the physical gauge boson masses. The details of the renormalisation procedure are given 
in section El 

Since the calculation involves the computation of more than thousand diagrams, it 
is convenient to employ computer-algebra tools. The generation of diagrams and Feyn- 
man amplitudes, including the counterterm contributions, was performed with the package 
FeynArts [25! ■ The program TwoCalc [2E! was applied for the algebraic evaluation of these 
amplitudes, which were reduced, by means of two-loop tensor-integral decompositions, to a 
set of standard scalar integrals. Throughout the calculation, a general gauge was used, 
and the gauge-parameter independence of the final result was checked algebraically. For 
the evaluation of the scalar one-loop integrals and the two-loop vacuum integrals we have 
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Figure 2: Two- loop vertex diagrams containing a triangle subgraph, which require a careful treat- 
ment of 75 in D dimensions. 

used analytical results as given in Refs. [27J|2E|, while the two- loop two-point integrals with 
non-vanishing external momentum have been evaluated numerically using one-dimensional 
integral representations with elementary functions j2H] • These allow a fast and stable calcu- 
lation of the integrals for general mass configurations. 

Since we use Dimensional Regularisation jHUJ|M] in our calculation, it is necessary to 
investigate the treatment of the Dirac algebra involving 75. It is known that a naively anti- 
commuting 75 respects all Ward identities of the Standard Model j22] ■ However, while it can 
safely be applied for all two- loop two-point contributions (for a discussion, see e.g. Ref. [T5] ) 
and most of the two-loop vertex- and box-type diagrams, it would yield an incorrect result 
for vertex diagrams containing a triangle subgraph (see Fig. EJ). This originates from an 
inconsistent treatment of the trace of 75 together with four Dirac matrices, which in four 
dimensions is given by Tr {757 /i 7 i/ 7 p 7 cr } = 4ie / " /p<J , while this trace would vanish when using 
the naively anti-commuting 75 in D dimensions. 

A mathematically consistent definition of 75 in D dimensions |31ll33] would require the 
introduction of additional counterterms to restore the Ward identities, which is a very tedious 
procedure at the two- loop level. For recent discussions on this topic, see Refs. [33J|33]. 

In order to calculate the class of diagrams in Fig. |2l we have first evaluated the triangle 
subgraph with a consistent 75 according to Refs. (here we made use of the package 

Tracer jSH] for checking). After adding appropriate counterterms to restore the Ward 
identities, the result differs from the result obtained using a naively anti-commuting 75 
only in terms proportional to the totally antisymmetric tensor & vpa , which are finite for 
D — > 4. Inserting this difference term into the two-loop diagrams, it turns out that the 
second loop only yields a finite contribution, so that it can be evaluated in four dimensions 
without further complications. After contraction with the external fermion line in the vertex 
diagrams, a non-zero contribution to the result for Ar is obtained from this term. 

3 On-shell Renormalisation 
3.1 One- loop renormalisation 

In the on-shell renormalisation scheme the mass parameters and coupling constants are 
related to physical observables. The gauge bosons of the U(l) and SU(2)l group, B^, W^' 2,3 
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are conveniently expressed in terms of their mass eigenstates, 




where W^, Z^ denote the fields of the massive vector bosons W, Z with masses Myv, Mz, A^ 
represents the massless photon field, and the weak mixing angle enters in the combination 

c w = cos # w = M w /M z , s w = sin 6 W = yjl - M^/M%. (6) 

The eigenstates of the Higgs doublet are the physical Higgs field H with mass Mh and 
the neutral and charged Goldstone bosons \i 'P^- We neglect mixing between the fermion 
generations throughout. 

In the following the conventions of [37| are adopted. In our approach all physical fields 
and masses as well as the electromagnetic coupling e are renormalised: 

Mi = Mi + 5Ml 

m fo = m f + 5m f, 
Z = {Z zz ) l l 2 Z + \5Z^A, (7) 

A = \bZ l7j Z + (Z 77 ) 1/2 A, 

/o L = (^ /L ) 1/2 / L (f = e,u e ,...), 
/o R = (z / R) i/ 2/ r 

Here the index indicates the bare quantities. The renormalisation constants for the masses, 
5Mx,Smf, the fields, Z x = 1 + 5Z X , and the charge, Z e = 1 + SZ e , are fixed by on-shell 
renormalisation conditions. 

The physical squared masses M\ are defined as the real part of the poles M.\ of the prop- 
agators D x . The field renormalisation constants are determined by demanding unity residues 
of the poles. This ensures that all Green's functions are finite and external wave function 
corrections need not be taken into account. With the index T denoting the transverse part 



e 



Z e e, 



M w + <5M W , 
Ml + 5 Ml , 

(Z WW2 w ± 



H n = (Z K W 2 H, 
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of the vector boson propagators, the on-shell renormalisation conditions read 
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The propagators are related to the one-particle-irreducible two-point functions T by 



D T W 



(Ft 

D u = -(f H )~\ 



( 



7Z 



D 
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pZZ fVyZ 
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\j + f f R + m f f f s ) \ 



(9) 
(10) 



-(f 7 ) 1 - - lF ^_i L 

Here f{ RS represent the left-/right-handed and scalar component of the fermion two-point 
functions, respectively. The two-point functions can be separated into a Born contribution 
and the self-energies, 

' * ' ' (11) 

(12) 



T^(k 2 ' 



(k 2 - M 2 )6 ab + Zf{k 2 ' 



i + e£, r (p* 



The hat indicates renormalised quantities, i. e. £ = E + counterterms. 

In addition to the aforementioned renormalisation conditions, one has the freedom to 
renormalise the Higgs tadpole t. Here the condition 



t + St = 



(13) 



is applied, which requires that all tadpole contributions are exactly cancelled by the coun- 
terterms 5t, so that no tadpoles need to be taken into account in the actual calculation. 
Using the renormalisation conditions eq. (jSJ) one obtains for the renormalisation constants 
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in eq. (J7J) at the one-loop level 



6M 2 , = Re{£W(M^)}, 
SM 2 = Re{£f (M z 2 )}, 



5M 2 =Re{E H (M 2 )}, 



5Z W 
5Z ZZ 

5Z^ Z 

5Z Z ~> 

SZ-rv 
5Z H 



Re{Ef(M w )}, 
Re{E| z/ (M z 2 )}, 



2 ^{^k^}' 



Ml 



S^ z (0) 
Ml ' 
-£^(0), 

-Re{£ H/ (M 2 )}, 



(14) 



5m y 



2 



m 



+ 2S S 
+ 2S 



(15) 



[m. 



Re{S{(m^)} -m) Rt 
= - Re{S^(m^)} - m) Re 

with Y, 1 {k 2 ) indicating the derivative of the self-energy with respect to k 2 . 

For the charge renormalisation an additional condition is required. Usually it is fixed 
by demanding that the electric charge e coincides with the coupling of the electromagnetic 
vertex in the Thomson limit, 



ieu(p)^^u{p). 



Employing the U(l) Ward identity this yields at one-loop order 



SZ K 



1 



5Z^ - 
2 2c w 



5Z Z \ 



(16) 



(17) 



The weak mixing angle is a derived quantity, expressed in terms of the gauge boson masses, 
see eq. (JBJ). Thus the renormalisation of M W ,M Z also determines the counterterm 5sw f° r 
sw 4j. At one-loop order one obtains 



<5s\y 



r 2 



2s 2 



w 



ml 



SM? 



w 



From two-loop order on, a one-loop sub-renormalisation is necessary for the Faddeev- 
Popov ghost sector, which is associated with the gauge-fixing part. It is possible to keep the 
gauge-fixing part invariant under renormalisation. For technical convenience, we arrange for 
this by a renormalisation of the gauge parameters in such a way that it precisely cancels the 
renormalisation of the parameters and fields in the gauge- fixing Lagrangian. 1 To this end, 

another approach to avoid counterterms emerging from the gauge-fixing sector, the gauge-fixing part 
could alternatively be added to the Lagrangian only after renormalisation. In this case the counterterms 
to the ghost sector arise only from the renormalisation of the gauge transformations and not from the 
gauge-fixing functions F v . 
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we start with the following, rather general form of the bare gauge-fixing term: 

£gf = -^((F"0 2 + (F z ) 2 + F + F- + F-F + ), (19) 

= (ayh^ + ^d^, (20) 

F z = (iiyh,Z» + ^d,A»-(gf M zX , (21) 

F ± = (^d^^i^M^^ (22) 

allowing two different bare gauge parameters for both W and Z, ^ 1 V,Z and £^' Z , and also 
mixing gauge parameters, £ 7Z and £ Z7 . The renormalised parameters shall comply with the 
i?£ gauge, providing one free gauge parameter for each gauge boson, £ 7 ,£ z ,£ w . With the 
following renormalisation prescription 




(23) 



' 2 M| 



Z* (24) 



~~ ^ _ ^ 2 ^ ' 

no counterterm contributions arise from the gauge-fixing sector. Here we have allowed for 
field renormalisation constants 5Z X , 5Z^ of the unphysical scalars \i ± as well. Starting at 
the two-loop level, counterterm contributions from the ghost sector have to be taken into 
account in the calculation of physical amplitudes. They follow from the variation of the 
gauge-fixing terms F a under infinitesimal gauge transformations, 56 b , b = 7, Z, ±, 

^ - a ^/f> ^/A&-W^'W. (26) 

These contributions can be derived from the action of the gauge transformations on the 
gauge and Goldstone fields as follows, 

a, - A li + d li ser' + ie(w+se--w-se + ), 

Z, - Z, + d,56 z -ie C ^(w+56--W;5e + ), 

> Wt + =f re (w£6ff> - ^W±56 z - A fJj 89 ± + ^Z,56±) , (27) 

» x -(Mz + -^^)^ z + -^(0 + (5r + 0-^ + ), 
V 2s w c w / 2s w v 1 



± 



X - 
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We have derived all the counterterms arising from the ghost sector (extending the results 
of Refs. |3*%l l3*9"] to a general gauge) and implemented them into the program FeynArts. 
In this way we could verify the finiteness of individual (gauge-parameter-dependent) build- 
ing blocks (e.g. the W- and the Z-boson self-energy) as a further check of the calculation. 
The explicit Feynman rules of the ghost sector including counterterms can be found in the 
appendix. 



3.2 Two-loop counterterms 



In the 0(a 2 ) calculation of the muon decay, two-loop counterterms arise for the transverse 
W propagator and the charged current vertex: 



w w 

*a/v\a,(vw"(><}a/vvvvw* 



SZ^ 2 - M w ) - 8Ml m - 5Z^ 1) 5M 2 Vt 



w. 



(i)' 



(2f 




y/2. 



1^- 



5s w 



+ - (8Z$ + 6Z™ + 5Z t 



(29) 



+ (1-loop renormalisation constants) 



The numbers in parentheses indicate the loop order. Throughout this paper, the two-loop 
contributions always include the subloop renormalisation. 

Concerning the mass renormalisation of unstable particles, from two-loop order on it 
makes a difference whether the mass is defined according to the real part of the complex 
pole of the S matrix, 



M = M — iMT, 



(30) 



or according to the pole of the real part of the propagator. In eq. (|5U|l M. denotes the 
complex pole of the S matrix as specified by the renormalisation conditions in eq. (jHJ. M, T 
are then interpreted as the corresponding mass and width of the unstable particle. For the 
real pole, on the other hand, we use the symbol M. It is determined by 



Re{(D T y\M 2 )} = 0. 



(31) 



In the context of the present calculation, these considerations are relevant for the renor- 
malisation of the gauge-boson masses, Mw and Mz- The two- loop mass counterterms ac- 
cording to the definition of the mass as the real part of the complex pole are obtained from 
eq. ©, 



5M 



W,(2) 



Re{s^ (2) (M w )} -5Ml Al) 5Zf x) +Im{s^ 1) (M w )}lm{sW (1) (M w )},(32) 



5M 



Z,(2) 



Re {E- )( M1)} - 5Ml A1) 5Z% + M ^ + M^*®) 



Ml 



+ Im {El^Ml)} Im {E^ 1} (M Z 2 )} . 



(33) 
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When compared with the mass counterterms according to the real-pole definition, SM^ 



(2) 



and <5M| there remains a finite difference, 



5M 
5M 



2 

W,(2) 
Z,(2) 



5M W,(2) 

5M Z 2 (2) 



Im{s^ 1) (M^)}lm{s?; (1) (M^)} 
+ Im {^(M 2 )} Im {S^ 1} (M Z 2 )} . 



(34) 
(35) 



It can easily be checked by direct computation that the difference terms in eqs. (jHlJ), (|33j) are 
gauge-parameter-dependent, thus showing that at least one of the two prescriptions leads to 
a gauge- dependent mass definition. The problem of a proper definition of unstable particles 
in gauge theories has already been addressed several times in the literature jlU] • However, the 
present work, for the first time, involves an explicit calculation of a physical process which is 
sensitive to the gauge-parameter dependent difference between the two mass renormalisation 
methods. In the previous results for Mw, incorporating terms up to 0(G 2 ^m 2 M%) [20j and 

Mn-dependent fermionic terms j2Ij, the contribution Im {e^ ^(M 2 ) j Im |St,(i)(M 2 )| was 
zero, making thus a strict distinction between the two mass definitions unnecessary at the 
considered order. 

Using a general gauge, we can test the two mass renormalisation prescriptions in 
our result by regarding the two-loop counterterms to physical observables, which should 
be gauge-parameter independent. In particular, we only find an invariant result for the 
counterterm to the weak mixing angle, <5s\y,(2), with the definition of the gauge-boson masses 
according to the complex pole. 

In order to verify the gauge-parameter independence of the mass counterterms 5M^ ^ , 
5M| / 2 n one needs an appropriate treatment of the Higgs tadpole diagrams, which do not 
contribute to physical observables. Alternatively to eq. (fT3*|) . one can include all Higgs tad- 
pole diagrams in the calculation by demanding the tadpole counterterm to be zero, 5t = 0. 
Technically, this corresponds to the inclusion of all tadpole diagrams not only in the two-loop 
self-energies, but also in the subloop renormalisation. In this case, also the mass countert- 
erms themselves are gauge-parameter independent when using the mass definition via the 
complex pole. 

These results confirm the expectations from S-matrix theory that the complex pole is a 
gauge-invariant quantity [4~Uj . 

We have thus adopted the complex-pole definition as given in eq. (J3~2j) and eq. (|3l!j) . Using 
this mass definition and expanding the gauge boson propagator around its pole 



(D T )-\ q 2 ) = (D T )-\M 2 ) + (q 2 -M 



_d_ 



(D T y\k 2 



k 2 =M 2 



M 



2\2 



one obtains with the renormalisation conditions eq. 



D T (q 2 



-i const. 



q 2 — M +iMT 



+ non-resonant terms, 



(36) 



(37) 



which corresponds to a Breit-Wigner parametrisation of the resonance line shape with a 
constant decay width. 
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Experimentally the gauge-boson masses are determined using a Breit-Wigner function 
with a running (energy-dependent) width, 



D ^ q2) K <r-AP +iq *T/M - (38) 

As a consequence of the different Breit-Wigner parametrisations, there is a numerical dif- 
ference between the experimental mass parameters (denoted as Mw, M z henceforth) and 
the mass parameters in our calculation, Mw, M z . The shift between these parameters is 
given by |12] M w ,z = M w ,z + r w z / (2M w ,z)- Since M w and M z enter on a different footing 
in our computation — M z is an experimental input parameter, while Mw is calculated - 
in order to evaluate the mass shifts we use the experimental value for the Z-boson width, 
T z = 2.944 ± 0.0024 GeV j^j, and the theoretical value for the W-boson width, which is 
given by Tw = 3G M M w /(2v2V)(l + 2a s /(37r)) in sufficiently good approximation. This 
results in M z « M z + 34.1 MeV and in the mass shifts M w ~ M w + 27.4 MeV and 
M w « M w + 27.0 MeV for M w = 80.4 GeV and M w = 80.2 GeV, respectively. 2 

For an extension of the renormalisation formalism for unstable particles to higher loop 
orders and to the field renormalisation of unstable particles, see Refs. |43j . However, in a 
physical process with particles in the initial and final state whose mass can be neglected, a 
treatment of complex poles is only necessary for internal particles. Since the field renormal- 
isation of internal particles does not contribute to the physical result, it is not necessary to 
examine this issue for our purposes. 

The two-loop charge renormalisation constant follows from the condition eq. (|16|). With 
the help of the U(l) Ward identity the electromagnetic vertex can be related to photonic 
two-point functions, thus resulting in the following relation between the renormalisation 
constants, which is valid in all orders of perturbation theory [39J: 

Z e ( (Z^ 2 + Sw + 5 f W S -^] = 1. (39) 
\ c w + £>c w 2 J 

Expansion up to 0(a 2 ) yields 

5Z e{2) = -\8Z^ - ^-5Z^ + (6Z e{1) ) 2 + I (5Z^ 2 - J-SZffi^y (40) 

Furthermore the two-loop field renormalisation constants for the external leptons are 
needed in eq. (I29|) . These can be easily obtained in the limit of vanishing masses and 
momenta of the external fermions, yielding 

*Zg = -Egj(mJ = 0). (41) 



2 The difference in T\y according to the way it is calculated, through the tree-level result parametrised with 
a, or the improved Born result parametrised with G^, or the improved Born result including QCD corrections 
(which is the one we used), is formally of higher order (i.e. beyond 0(a 2 )) in the calculation of Mw- Its 
numerical effect is nevertheless not completely negligible; it changes the shift in Mw by about —2.9 MeV if 
the tree- level result for Tw parametrised with a is used and by about —1.4 MeV if the G M parametrisation 
of the Born width (without QCD corrections) is employed. 
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Figure 3: Virtual (a) and real (b) QED corrections to muon decay in the Fermi Model. 

4 Extraction of Fermi Model QED contributions 

In the evaluation of Ar, the IR-divergent QED corrections that are already contained in the 
Fermi Model QED factor have to be extracted. For the two-loop calculation presented here, 
the corresponding Fermi Model contributions consist of virtual and real photonic corrections 
of order 0{ot) and of order 0(a 2 ) with one closed fermion loop, see Fig. where it is 
understood that all lepton and quark flavours can appear in the loop. Denoting the virtual 
corrections to the Fermi Model by Agv and the real corrections by AgR, this reads 

l-MFcrmil 2 = \M Born \ 2 (1 + Aq) = \M Bo m\ 2 (l + Ag^ + Aqg> + Ag^ + Ag^ 2) ) . (42) 

The calculation of the virtual corrections to muon decay in the full Standard Model 
involves box-type diagrams with IR divergences. In the following, all Standard Model con- 
tributions involving photons in the loop are encompassed by the quantity Ar. The one-loop 
QED corrections At v originate from the diagram given in Fig. |H(a). At two-loop order one 
can distinguish between corrections with only electromagnetic couplings in addition to the 

( 2 \ 

tree-level couplings, Ary^, see Fig. IH(b), and corrections with additional QED and non- 

QED couplings, Ary e ^ weak , see Fig. |U(c). Here it is always understood that the two-loop 
corrections involve one closed fermion loop. By performing a tensor integral decomposition 
of these classes of diagrams one observes that they can be expressed in terms of the virtual 
corrections in the Fermi Model and an IR-fmite remainder Arf r , 

Ag(r } +2Ar{ r a) (43) 
Aqf ) + 2Arg ) (44) 
2 (Ag^ + Ari Q) ) Ar£L + 2 Ar£j. (45) 

Note that the factor 2 in these formulae arises due to fact that Ag enters linearly into the 
muon decay width, see eq. (JTJ), while there is a quadratic dependence on Ar, see eq. (jSJ). The 



Ar^ 
Ar (a2) 
Ar (a2) 

V, em/weak 
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Figure 4: Examples for virtual IR-divergent diagrams contributing to muon decay in the Standard 
Model. Besides the one-loop diagram (a), at two- loop order, there are diagrams with four (b) and 
two (c) electromagnetic couplings. 




Figure 5: Examples for real bremsstrahlung diagrams in the two-loop calculation of muon decay 
in the Standard Model, involving one-loop diagrams (a) and two-loop diagrams with four (b) and 
two (c) QED couplings. 

finite remainders are then combined with all remaining virtual Standard Model contributions 
into Ar^ a \ Ar^ a2 \ In eq. (I45J1 Ar^ m corresponds to the non-QED one-loop corrections with 
a closed fermion loop. 

Besides box-type diagrams, IR divergences are also present in the field renormalisation 
of the external leptons. Here the correspondence to the Fermi Model contributions is trivial. 

Similar to the virtual diagrams, the real bremsstrahlung corrections to muon decay in 
the Standard Model can be divided into the one-loop contribution Tp, , two-loop corrections 

( 2 \ 

with additional electromagnetic couplings only, r^^, and mixed QED/non-QED couplings, 

( 2 \ 

r Rem/weak' see Fig- El Exploiting the fact that the momentum transfer q 2 through the W 

boson propagator of these diagrams is much smaller than the W mass, q 2 <C M|r, in the 
limit of zero momentum transfer the real contributions can be reduced to the real Fermi 
Model contributions, 

Ar£" = Aq^ (46) 

Ar£2 = Agf > (47) 
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Figure 6: In order to extract the QED corrections already present in the Fermi Model from the 
Standard Model computation of Ar, differences between QED loop diagrams in the Standard Model 
and Fermi Model of the same order have to be evaluated, here exemplified for the one-loop case. 



At (q2) - 2 A« (a) Ar (a) (At 

^ T R,em/weak — Z ^R ^'ferm- 

In total, the contributions to the two-loop Standard Model matrix element amount to 



|.Msm| 2 = l-MBornl 2 



1 + Ar<«> + Ar(° 2 )) 2 + AgW + Aqf ] + 2 Aq^ Ar^ m 

+ AgW + Agf)+2Ag^Ar{ e i 

3 



(49) 



^Bo m | 2 [(l + Ag)(l + Ar) 2 + 0(« 



which can be written in the factorised form at least up to two-loop order including one closed 
fermion loop. Contributions with two closed fermion loops at 0(a 2 ) are not present in the 
Fermi Model and do not contain any IR divergences. Comparing eq. (jlSj) with eq. (}42|) one 
obtains 

|A^SM| 2 = lA^Fcrmil 2 (1 + Ar) 2 , (50) 

showing that a factorisation of electromagnetic corrections to the Fermi Model and the 
remaining electroweak corrections in the Standard Model according to eq. (JIJ, (J2J) is possible 
at least up to the given order. 

The calculation of the remaining terms Arf r , in which the QED Standard Model and 
Fermi Model contributions in eq. I)4H |) -(|45 J) differ, requires the subtraction of Fermi Model 
diagrams from Standard Model graphs, as shown in Fig. El These terms are IR-finite but UV- 
divergent and therefore require regular isat ion. In our approach, dimensional regularisation 
turns out to be problematic for this purpose since for the computation of the fermion lines 
in diagrams like those in Fig. El we use the Chisholm identity 

IpTlvlp = -K^ P *l CT l<o + Qnvlp + 9u P lp ~ 9pplv (51) 

This identity, however, is only valid in 4 dimensions. In order to circumvent this problem 
we employ Pauli-Villars regularisation for the QED corrections to the Fermi vertex. The 
combination of these vertex corrections, Fig. El (a), with the QED part of the field renor- 
malisation of the external leptons forms an UV-fmite quantity. It is therefore possible to 
evaluate this combination using Pauli-Villars regularisation (PaVi) and employ dimensional 
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regularisation (DReg) for the rest. This can formally be written as 



Arf r = [QED box graphs in SM] 
+ 



DReg 



[QED vertex corr. to FM] 



PaVi 



DReg 



2^em + 2^em 



(52) 



PaVi 



The index "em" at the field renormalisation constants SZjt^, SZ^ indicates that only the 
QED-like diagrams of the lepton self-energies are taken into account for the calculation of 
these constants. Since the Standard Model box diagrams are UV finite, they can also be 
computed with a Pauli-Villars regulator. Thus the cancellation of the IR divergences between 
the two terms in the first line of eq. (p)2"j) proceeds in a straightforward manner and no IR 
regulator is required. 

A similar cancellation of IR divergences takes place for the terms in the second line of 
eq. (J52j) . This can be made explicit by introducing the Pauli-Villars regulator A in photon 
propagators according to the replacement 



1 



1 



1 



k 2 k 2 - A 2 



(53) 



with k being the photon momentum. In the difference in the second line of eq. (p)2"j) this 
corresponds to the replacement of the photon propagator in the lepton self-energies by 



1 

k 2 



1 

Y 2 



l 

k 2 



k 2 - A 2 



k 2 - k 2 '' 



(54) 



or, for the case of two photon propagators in the loop, 

1 



1 1 

k 2 X k 2 



1 1 

fc2 X W ~ 
1 

9 x 

k 2 k 2 



1 

k 2 



k 2 -A 2 



x 



1 

V 2 
1 



k 2 - A 2 



A s 



x 



k 2 - A 2 k 2 - A 2 



(55) 



It can be seen that this replacement effectively leads to the introduction of massive photons 
with mass A so that no IR divergences are present anymore. 



5 Numerical Results 

We shall now discuss the numerical evaluation of our result for Ar. It should be noted 
that our definition of Ar according to eq. (0) is based on the expanded form (1 + Ar) 
with Ar = Ar^ + Ar^ Q - 1 + . . . rather than on the resummed form 1/(1 — Ar). The terms 
obtained at two-loop order from a resummation of leading one-loop contributions are directly 
contained in our two-loop contribution to Ar. The following contributions to Ar are taken 
into account 

^ r _ A r («) + /\A aa ^ + Ar^""^ + Ar^'" 2 ' + Ar'^" 2 ', (56) 

where Ar'"' is the one-loop result, eq. ©, Ar*"" 5 ) and Ar^ aa ^ are the two-loop [H] and 
three-loop QCD corrections, while Ar^ Nta ' is the new electroweak two-loop result. 
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Figure 7: Two-loop contributions with one closed fermion loop to Ar as a function of the Higgs 
mass. The plot shows the full result, Ar^ Nfa \ the contributions from the top-bottom doublet, 
Ar^*" >, and the light-fermion doublets, Ar^ Nlia \ as well as the light-fermionic contribution 
with subtracted resummed terms proportional to Act, Ar^ f "_] oop , see eq. (|58|). 

The symbolic notation (N^a 2 ) encompasses the contribution of all diagrams containing one 
fermion loop, i. e. both the top/bottom and light-fermion contributions. Correspondingly, 
the term Ar^f 2 ™ 2 ' contains the pure fermion-loop contributions in two-loop order. 

The pure fermion-loop contributions in three- and four-loop order turn out to be numer- 
ically small, as a consequence of accidental numerical cancellations, with a net effect of only 
about 1 MeV in M w (using the real-pole definition of the gauge-boson masses) [T2J. Fur- 
thermore, also the leading three-loop contributions for a large top-quark mass in the limit of 
zero Higgs mass, proportional to a 3 m^ and a 2 a s mf, have very little impact on the prediction 
of M\y [13J. Therefore these two corrections have not been included in this analysis. 

In Figure [7| and Table ^ numerical values for different two- loop contributions with one 
closed fermion are given as a function of the Higgs boson mass Mh, using M\y = 80.451 jB], 
m t = 174.3 03] and Aa = 0.05911 j3H]- The contributions with two closed fermion loops 
are not given in the figure and table since they are independent of Mh. Numerically, they 
yield a contribution of Ar' w f ° ' = 16.3 ■ 10~ 4 . The Higgs-mass dependence of the two-loop 
result for Ar agrees perfectly with the result previously obtained in Ref. [23]. 

It can be seen that both corrections with a top-/bottom-ioop, Ar^" ', and with a 
light-fermion loop, Ar^ 1 '" ' yield important contributions. At first glance it looks surprising 
that the light-fermion contributions even dominate over the top-/bottom-contributions for 
large Higgs masses (Mh ^ 300 GeV), which seems to endanger the validity of the large-mt 
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M H / GeV 


Ar (AT f a 2 )/ 10 -4 


Ar (N th a 2 ) /10" 4 


Ar-Wf" 2 ) - 2AaAr£> /10~ 4 


65 


20.1 


15.7 


2.0 


100 


20.8 


14.7 


2.2 


300 


23.3 


11.7 


3.1 


600 


24.9 


8.9 


3.9 


1000 


24.5 


5.1 


4.6 



Table 1: Two-loop contributions with one closed fermion loop to Ar for different values 
of the Higgs mass. Besides the full result, Ar^ fa \ also the contributions from the top- 
bottom doublet, Ar( Ntha and the light-fermion doublets, Ar^ lfQ \ are given. From the 
latter the term proportional to Aa originating from the resummation prescription eq. (|57j) 
is subtracted. 



expansions (TBIEIIIISI]- However, in previous analyses, resummation prescriptions ^T] have 
been derived in order to obtain partial terms of the two- loop result. With the replacement 



1 + Ar 



1 



1 - Ar 



(57) 



the term 2AaAr[^, generated from the charge renormalisation in bosonic one-loop terms, 
is correctly predicted jTTj, as we have checked by comparing with the full result, Ar^ Nua 
Therefore, in order to demonstrate the effect of the new two-loop contribution, the difference 



sub, 1— loop 



Ar m^) _ 2A«Ar£ 



(5f 



is shown in Figure [3 and Tabled This expression does not exceed the top-/bottom con- 
tributions for any value of the Higgs mass below 1 TeV. The new contribution to Ar from 
diagrams with a light-fermion loop amounts up to 3.3 ■ 10~ 4 which corresponds to a shift in 
M w of > 5 MeV. 

The prediction for Mw is obtained from eq. (0) by means of an iterative procedure, since 
Ar itself depends on M w , 



i 



My < - + A 

^2 Y V2G,Mi 



Ar(M w ,M z ,M H ,m t ,...) 



(59) 



In Fig. |H]the prediction for Mw based on the results of eq. (JSTIj) is shown as a function of Mh 
for m t = 174.3 ±5.1 GeV jUj and Aa = 0.05911 ±0.00036 For comparison, the present 
experimental value, M^ p = 80.451 ± 0.033 GeV [Sj, and the experimental 95% C.L. lower 
bound on Mh (Mh = 114.1 GeV |1B]) from the direct search are also indicated. The plot 
exhibits the well-known preference for a light Higgs boson within the SM. In particular, the 
theoretical prediction (including the band from a variation of m t , which at present dominates 
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Figure 8: The SM prediction for J% as a function of Mr for m t = 174.3 ±5.1 GeV is compared 
with the current experimental value, M^ p = 80.451 db 0.033 GeV and the experimental 95% 
C.L. lower bound on the Higgs-boson mass, Mh = 114.1 GeV [io] . 



the uncertainty of the prediction, see below, and Aa within la) cannot be matched with the 
la region of M^ p and the 95% C.L. exclusion limit on Mh. 

We have compared our results with those of an expansion for asymptotically large values 
of m t up to (9(G^m^M|) [2UII1Z1- The results are shown in Tabled for different values of 
Mh. The values for the input parameters are taken from Ref. [201, i-e. m t = 175 GeV, Mi = 
91.1863 GeV, Aa = 0.0594, a s (Mz) = 0.118. One observes a relatively good agreement, 
with maximal deviations in Mw of about 5 MeV. 

It should be noted that the deviations in the last column of Table |2]can not be attributed 
solely to differences in the two-loop fermionic contributions, because the results also differ 
by a slightly different treatment of higher-order terms that are not yet under control. 

In a further analysis, we have aimed at reducing the latter deviations as far as possible 
in order to focus on the effects from the two-loop top-quark and light-fermion contributions 

(see also the discussion in Ref. |22J)- While the result of Ref. j2Dl contains a term (^Ar^j 
generated from the purely bosonic one-loop contributions by the resummation eq. 1)5 70. no 
such term is included in our result. A second possible source of deviation could be caused by 
different implementations of the QCD corrections. We have therefore performed a compari- 
son in which the QCD corrections were removed from both results. With these modifications 
the maximum deviation between the results is not decreased, see second column in Table HH 
but the maximal difference in the Higgs-mass dependence Mw(Mh) — Mw(Mh = 65 GeV) 
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M H / GeV 


M w / GeV 


M^ pa / GeV 


AM W / MeV 


65 


80.3985 


80.4039 


-5.4 


100 


80.3759 


80.3805 


-4.6 


300 


80.3039 


80.3061 


-2.2 


600 


80.2509 


80.2521 


-1.2 


1000 


80.2122 


80.2129 


-0.7 



Table 2: The two-loop result for Mw based on eq. ()56|) is compared with the results of an 
expansion in m t up to (9(G^m^M|) [2DEZ1, M^ pa . The last column indicates the difference 
between the two results. 



Ma/ GeV 


AM W / MeV 


AM W / MeV 


65 


-5.4 


-2.3 


100 


-4.9 


-1.6 


300 


-4.0 


0.0 


600 


-4.3 


0.5 


1000 


-5.0 


0.5 



Table 3: Investigation of different sources of deviations between the two-loop result for Mw 
based on eq. (jSUjl and the results of j20EZl, M w cpa . In the second column (AM W ), differences 
in the QCD implementation and in the resummation of bosonic one-loop terms have been 
eliminated. In the third column (AM W ), in addition the two-loop contributions from light 
fermions are excluded. 

is reduced from 4.7 MeV in Table H to 1.4 MeV. 

It is also interesting to separately examine the effects of the top-bottom contributions that 
are not contained in Ref. [2U] an d of the two-loop terms from the remaining light-fermionic 
flavours. In the third column in Table E3 we have therefore excluded all light-fermionic 
0(a 2 ) contributions from the comparison. This is achieved by subtracting the expression 
Ar^"" 2 ' — 2Ar| f °' ) Ar^g, where the second term 2Ar^ Ar^ accounts for the light-fermionic 
terms that were included in Ref. [20] by means of the resummation presciption eq. ([57)1 . The 
remaining deviations between the results, which now only contain top-/bottom contributions 
at the two-loop level, are somewhat smaller, while there are larger differences in the Higgs 
mass dependence of up to 2.8 MeV. 

In Ref. P a simple formula was given which parametrises our full result for Mw, 

M w = M w - a dH - c 5 dH 2 + c 6 dH 4 -c 2 da + c 3 dt - c 7 dH dt - c 4 da s , (60) 



19 



where 



dH = ln llOOGev) ' dt = {lUjGev) d ° = 0^924^ ^ = Tll^- 1 ' (61) 

and M z = 91.1875 GeV [3] has been used. By a least square fit we have obtained for the 
coefficients ci,...,Oj 

Myj- = 80.3755 GeV, c 4 = 0.0763 GeV, 
ci = 0.0561 GeV, c 5 = 0.00936 GeV, 

(62) 

c 2 = 1.081 GeV, c 6 = 0.000546 GeV, 

c 3 = 0.5235 GeV, c 7 = 0.00573 GeV. 

The parametrisation of eq. (ffiUj) approximates our full result for Mw within 0.4 MeV for 
65 GeV < Mh < 1 TeV and the other input values within their la experimental bounds. 

Since this region of validity is in general not sufficient for global fits of the Standard 
Model, here we supply a more elaborate parametrisation, including the dependence on the 
Z-boson mass, 

M w = M^-d 1 dR-d 2 dR 2 + d 3 dR 4 -d i da+d 5 dt-d 6 dt 2 -d 7 dRdt-d s da s + dgdZ, (63) 

with dZ = Mz/(91.1875 GeV) — 1 and the other symbols as given in eq. ([61)1 . With the 
following values for the coefficients di, . . . , dg, 





= 80.3756 GeV, 


dr. 


= 0.5236 GeV, 


d\ 


= 0.05619 GeV, 


d 6 


= 0.0727 GeV, 


d 2 


= 0.009305 GeV, 


d-7 


= 0.00544 GeV, 


d 3 


= 0.0005365 GeV, 


ds 


= 0.0765 GeV, 


d.i 


= 1.078 GeV, 


dg 


= 115.0 GeV, 



(64) 



the full result for Mw is approximated by eq. (J63|) better than 0.3 MeV for 65 GeV < Mh < 
1 TeV and 2a variations of all other experimental input values. 



6 Remaining theoretical uncertainties 

Presently, the prediction of the W mass from Ar is mainly affected by the experimental error 
in the top mass determination, m t = 174.3 ± 5.1 03]. This induces an error of ~ 30 MeV 
in the predicted W mass. It is expected that the LHC can reduce the error on the top mass 
down to about 1.5 GeV [3B1 and a high-luminosity linear collider even to below 200 MeV jH], 
resulting in an error in the Mw-prediction from the m t -uncertainty of ~ 10 MeV and ~ 1.2 
MeV, respectively. Another important source of uncertainty is the experimental error in the 
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determination of Aa, which in a recent analysis was quoted to be 36 • 10~ 5 inducing an 
error of ~ 6.5 MeV in the predicted W mass. It is expected that this uncertainty will be 
further reduced significantly in the future • On the other hand, the experimental error of 
the direct measurement of the W mass, currently 33 MeV, is expected to reduce to 15 MeV 
for the LHC [7j and 6 MeV for a linear collider running at the W pair threshold 

Concerning the theoretical prediction, there are three main sources for uncertainties 
induced by unknown higher orders: the missing purely bosonic two-loop contributions, three- 
loop electroweak contributions and the lowest missing QCD corrections of order 0(a 2 a s ) and 
0{aal). 

For the three-loop O (a 3 ) electroweak corrections, partial results are available. In par- 
ticular the contribution from purely fermionic loops J2| is known, which amounts to a net 
effect of about 1 MeV in M w . Recently, the leading terms for large top masses proportional 
to m® [13] . which enter via the quantity Ap, have been calculated, having an effect of much 
less than 1 MeV on Mw- However, for the case of the 0(a 2 ) correction it turned out that 
the formally leading term oc m\ JH] i n the limit Mh = is suppressed relative to the next 
term oc m\ [20J, which suggests that the full 0(a 3 ) corrections could be of (9(1 MeV). Al- 
ternatively, one could try to estimate their size from residual scheme dependencies of the 
0(a 2 ) result. For example, the W width is needed at tree-level in order to translate be- 
tween the different Breit-Wigner parametrisations mentioned in section l3~2l Whether T^y is 
parametrised with a or G M is formally of order 0(a 3 ) and also results in a shift in Mw of 
0(1 MeV). 

For the QCD correction of order 0(a 2 a s ), the leading contribution oc m\ in an expansion 
for large m t has been calculated in the limit of vanishing Higgs mass ^3J. It turned out to 
result in a W mass shift of only less than 0.5 MeV. However, as explained above, the formally 
leading term oc vn\ can be suppressed relative to the sub-leading terms, so that the total 
0(a 2 a s ) contribution could be considerably larger. 

Higher order QCD contributions may be estimated from the renormalisation scale de- 
pendence of the available results. For this purpose running MS values at different scales are 
used for the top-mass m t or strong coupling a s . With this variation in the 0(a 2 ) result we 
obtain a shift in Mw of m 3.8 MeV, which gives an estimate of the 0(a 2 a s ) contributions. 
Accordingly, from the scale dependence of the 0(aa 2 ) result, we estimate the effect of the 
0(ac%) term to yield « 0.7 MeV. 

A second method for estimating the missing QCD corrections relies on the assumption 
that the ratios of consecutive coefficients in the pertubative series do not change very much. 
In this case the ratio between the 0(a 2 a s ) and 0(a 2 ) contributions to Ar should be of the 
same size as the ratio between Ar( aas ) and Ar' a ' 3 . From this we deduce that the effect 
of the 0(a 2 a s ) contribution on Mw is about 3.5 MeV, which is in nice agreement with the 
previous estimate. In a similar manner, we can compare the ratios Ar^ aa - )/Ar ( - aa ^ and 
^Kl^j-K) aric [ obtain « 0.7 MeV for the impact of the 0(aa 3 ) contribution. 

The third important higher order contribution, the bosonic 0(a 2 ) correction, is more 
difficult to estimate. In a simple approach the one-loop bosonic contribution is resummed 

3 This ratio amounts to » 12% in accordance with a s (Mz). 
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Figure 9: Theoretical error band of the M\y-prediction due to unknown higher order contributions. 



according to the prescription eq. (|57jl which results in a term ^Ar^J . This term shifts the 
W mass by less than 0.5 MeV for M H = 100 GeV, but more than 2.5 MeV for M H = 1 TeV. 

We obtain the total theory uncertainty by linearly adding up all sources for theoretical 
errors. This results in an uncertainty for M\y of 6 MeV for light Higgs masses and about 8 
MeV for Mh ~ 1 TeV, which is similar to the value given in Ref. jSHj- The error band due 
to these theoretical uncertainties is shown in Figure El 

Recently the new two-loop results for the prediction of the W-boson mass have been 
implemented into the Standard Model fits with Zfitter and are subject to the latest 
LEP electroweak analyses [5]. While the effect on the predicted value for M w is relatively 
small compared to the experimental error, it induces a significant shift in the prediction of 
the effective leptonic weak mixing angle sin 2 9^ according to 



where k incorporates the contributions from radiative corrections. The effect of inserting 
the new result for M\v in eq. (J55j) instead of the previous result obtained from an expansion 
in powers of m t j2U| amounts to an upward shift of about 8 ■ 10 -5 , which is about half the 
experimental error of 17 ■ 10~ 5 [5J. Since the corresponding complete fermionic two-loop 
corrections for k are not yet known, this shift has been treated as a theoretical uncertainty 
and is represented as a rather wide band in the well-known blue-band plot jS]. 





(65) 
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Figure 10: Examples for types of bosonic two- loop diagrams with internal Higgs bosons contribut- 
ing to muon decay. 

7 Higgs-mass dependence of bosonic two-loop result 

As a first step towards a full 0(a 2 ) result for Ar we have calculated the dependence of the 
bosonic two-loop corrections on the Higgs-boson mass. This includes the evaluation of all 
diagrams without closed fermion loops which contain internal Higgs bosons or Mn-dependent 
scalar couplings. Some typical examples are given in Figure [OB This subset of the complete 
bosonic two-loop corrections can be evaluated with the methods described in sections 
In particular, the factorisation of IR-divergent QED corrections as in eq. (|49|) also applies 
for the bosonic Mn-dependent contributions. 

In order to study the Higgs-mass dependence, the subtracted quantity 

Arg su b(M H , M H ) = Ar£>(M u ) - Ar£P(M°) (66) 

is considered, using a fixed offset value for the Higgs-boson mass, M H = 100 GeV. The 
contribution of the Mn-dependent diagrams to Ar^" s ' su1) (Mh, M h ) forms a finite and gauge- 
parameter independent quantity, as we have explicitly checked. This analysis is in analogy 
to Ref. \21\ . were the corresponding quantity for the fermionic two-loop contributions was 
studied. 

In Table |3] the variation of the prediction for the W-mass as a function of the Higgs 
mass Mh is shown without and with the bosonic two-loop terms, using the input values of 
Table ^ As before the values are given in terms of the subtracted quantity 

M w , S ub(M H ) = M W (M H ) - M W (M H = 100 GeV). (67) 

Figure ITT1 shows how the slope of the Higgs-mass dependence is modified due to the inclusion 
of the bosonic two-loop two-loop terms. The maximum change amounts to less than 2 MeV 
in the region 100 GeV < Mh < 1 TeV. As a consequence, from the Mn-dependence we get 
no indications for any particularly large effects in the full bosonic two-loop corrections to Ar. 
In this context we would like to point out the observation that the Higgs-mass dependence 
of the fermionic two-loop corrections [21] provides a rough assessment of the effect of the full 
two-loop corrections [1. This supports the estimation in section |U] that the expected size of 
the purely bosonic 0(a 2 ) contributions is relatively small. 
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M H / GeV 


M#7 ub /MeV 


M^ u + bos /MeV 


AM w , S ub/MeV 


100 











200 


-43.1 


-42.6 


0.5 


400 


-93.7 


-93.0 


0.8 


600 


-124.7 


-123.8 


1.0 


1000 


-163.4 


-161.7 


1.7 



Table 4: Shift in the predicted W mass caused by varying Higgs-boson mass M H when 
including fermionic two-loop corrections (Mw™ b ) and all two-loop corrections ^M-^ r ™^ bos ^ . 
AMw,sub gives the difference between these two contributions, i. e. the effect of the bosonic 
loops. 



1.75 




200 400 600 800 1000 

M H [GeV] 



Figure 11: Variation of the Higgs-mass dependence of the Mw-prediction due to the inclusion of 
bosonic two-loop corrections. 
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Figure 12: Comparison of prediction and measurement for Mw using expected experimental un- 
certainties at the LHC and the current central values. 



8 Prospects for future colliders 

In the following we illustrate the accuracy that can be reached with future colliders concern- 
ing tests of electroweak physics. Taking the current central values of the experimental input 
values, Fig. II 21 shows the situation that can be obtained with the LHC with expected errors 
for Mw and m t of 8Myj = 15 MeV and 5m t = 1.5 GeV, respectively. Even more impressive 
results could be achieved by a high-luminosity linear collider running at low energies, where 
errors of SMy? = 6 MeV and 5m t = 200 MeV may be obtained jH], see Fig. ED in both fig- 
ures we furthermore assumed that the error in the shift of the electromagnetic fine structure 
constant, Aa, will be cut to half of the present value. 



9 Conclusion 

In this paper, the evaluation of the complete fermionic two-loop contributions to the M^- 
Mz mass correlation was described, elucidating the applied techniques and the implications 
of the new result. 

The renormalisation within the on-shell scheme was described in detail. In particular, 
the definition of the gauge-boson masses via the complex pole of the S matrix was studied, 
ensuring in particular gauge-parameter independence of the renormalised weak mixing angle 
and the gauge boson masses. The latter requires to take tadpole contributions into account. 
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Figure 13: Comparison of prediction and measurement for M\y for expected experimental errors 
at a high-luminosity linear collider, assuming current central values. 



It was shown how the radiative corrections in the Standard Model can be factorised from 
the QED corrections within the Fermi Model. 

The result for Mw was expressed in terms of an accurate numerical parametrisation valid 
for all values of the Higgs mass up to 1 TeV. A detailed comparison with a previous result 
obtained by an expansion in powers of m t up to next-to-leading order was performed. Here 
the effects of the top/bottom and light-fermion contributions were studied separately and 
found to yield a contribution of a few MeV to the prediction of M w each. 

Furthermore, the remaining theoretical uncertainties due to unknown higher orders were 
discussed and an overall uncertainty of the W-boson mass prediction of ~ 6 MeV was 
estimated for light Higgs-boson masses. A careful treatment of the theoretical uncertainties 
proved to be important for precision tests of the Standard Model. The situation for present 
experimental uncertainties was contrasted to the capabilities of aspired future colliders. 

As an additional result, the Higgs-mass dependence of the purely bosonic electroweak 
two-loop contributions was computed, so that the only yet missing piece of the complete 
two-loop calculation for Ar, i. e. muon decay, is a constant (Mn-independent) contribution. 
The numerical impact of the bosonic two-loop corrections on the Higgs-mass dependence 
of the M w -prediction is relatively small, in accordance with our estimates for theoretical 
uncertainties. 
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Appendix 

In the appendix the explicit Feynman rules for the ghost interactions in the Standard Model 
are listed. In the vertices all particles are considered as incoming. The following expressions 
comply with a non-renormalisation of the gauge-fixing sector and the use of a linear 
gauge, introducing the (renormalised) gauge parameters £ 7 ,£ z ,£ w according to eqs. (J33J)- 
(J25|) . Counterterms in the one-loop approximation are included. The Feynman rules for the 
other sectors of the theory can be found e. g. in Ref. |57| . 
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